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ABSTRACT 


In  this  report  first  are  summarized  the  works  carried 
out  on  unsteady  effects  in  droplet  evaporation  and  combus- 
tion and  on  the  behaviour  of  hydrazine  catalytic  thrusters 
working  under  steady  state  conditions.  The  research  performed 
in  these  two  fields  was  presented  in  more  detail  in  previo\ s 
scientific  reports. 

The  third  topic,  that  is  presented  in  more  detail, 
is  on  the  influence  of  the  temperatures  of  liquid  hydrazine 
and  catalyst  particle  temperatures  on  the  vacuum  start 
characteristics  of  a reactor.  The  model  used  assumes  the 
existence  of  a liquid  front  that  first  advances  towards 
the  interior  of  the  catalyst  particle  by  porous  action 
and  then  recedes  towards  the  exterior.  An  order  of  magni- 
tude estimate  seems  to  indicate  that  during  the  process, 
convective  heat  transfer  dominates  over  conduction,  and 
the  analysis  shows  the  existence  of  a thin  thermal  front, 
whose  structure  is  calculated.  Then,  the  liquid  penetration 
is  only  governed  by  the  initial  particle  temperature.  For 
low  enough  initial  particle  temperatures  the  thermal  front 
structure  reaches  the  liquid  front  and  the  liquid  penetra- 
tion is  also  influenced  by  the  initial  hydrazine  temperature. 


ftSKFlSAllWf.,., 


Ui"  ni 


. BftTBISilTlOMyKUSlilp  j 
Sitt.  < if' ,-v 


' i i 




1 


1.  IKTR0DUCC10N  AND  SUMMARY  OF  PREVIOUS  v'?, SEARCH  WORK 


During  the  report  period  work  has  been  done  on 
the  f jllowing  three  areas. 

a)  Unsteady  Effects  in  Droplet  Evaporation  and 
Combustion. 

b)  Steady  State  Analysis  of  Hydrazine  Catalytic 
Thrusters . 

c)  Influence  of  both  Liquid  Hydrazine  and  Catalyst 
Particle  Temperatures  on  the  Vacuum  Start  Characteristics 
of  a Reactor. 

The  work  carried  out  on  the  first  two  subjects 
has  already  been  presented  in  detail  in  three  previous 
Scientific  Reports  and  only  will  be  summarized  in  the 
present  report.  Work  on  topic  c)  has  been  done  mainly 
during  this  year  and  will  be  presented  in  more  decail. 


1.1  Unsteady  Effects  in  Droplet  Evaporation  and  Combustion 


The  work  carried  out  on  this  subject  has  been 
presented  already  in  a report  by  Crespo  and  Linan  this 
has  also  appeared  with  some  modifications  as  a journal 
article  (see  a.Lso  ref.l). 

The  problem  considered  was  that  of  the  droplet 
vaporization  and  combustion  in  a stagnant  atmosphere 
taking  into  account  the  unsteady  effects  in  the  gas 
phase.  Previous  works  on  the  subject  assumed  the  process 
to  be  quasisteady,  based  on  the  fact  that  the  ratio 


of  the  densities  of  the  ambient  gas  and  of  the  liquid 
was  small.  In  this  work  a perturbation  analysis  is  made 
in  which,  by  considering  the  above  mentioned  ratio  as 
small,  the  asymptotic  nature  of  the  quasisteady  (QS) 
solutions  is  shown.  The  analysis  shows  that  corrections 
to  the  QS  theory  are  of  the  order  of  the  square  root 
of  the  ratio  of  the  densities  and  that  the  solution 
is  non-steady  (NS)  at  distances  from  the  droplet  of 
the  order  of  the  droplet  radius  divided  by  the  square 
root  of  the  ratio  of  the  densities. 

It  is  found  that  the  droplet  lifetime  is  given 
by  the  simple  analytical  expression 


!em0/ir 


Where  t^  is  the  lifetime  that  would  be  obtained  from 

the  QS  analysis,  e is  the  abovementioned  ratio  of  the 

densities  of  gas  and  liquid  fuel,  and  m^  is  a parameter 

defined  in  ref.  1 and  that  is  inversely  proportional 

co  t^Q.  A comparison  is  given  in  figure  1 of  the  result 

given  by  expression  (1)  with  those  obtained  numerically 

2 

by  Duda  and  Vrentas  for  the  problem  of  disolution  of 
spheres,  that  can  be  shown  to  be  mathematically  similar 
to  this  one.  The  agreement  is  good,  even  for  values 
of  en>Q  of  order  unity,  for  which  the  present  perturbation 
theory  can  not  be  applied. 


For  the  combustion  rase  the  flame  position  and 
extinctiot'  time  have  also  been  investigated.  When  the 
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stoichiometric  ratio  oxidizer  fuel  is  not  very  large 
it  is  found  that  QS  theories  give  flame  position 
correctly  with  a relative  error  of  order  /e\  In  many 


practical  cases  this  ratio  is  large,  even  of  the  order 
-1  /? 

~ C ~ ' ''  J Cl - * ~ *1 i J a.1 J „ 1 >tr»  • 


of  e ",  and  the  flatre  is  located  then  in  the  HS  region 
and  its  position  is  not  given,  even  in  first  approximation, 


by  the  QS  theory.  Comparisons  have  been  carried  out 

* A.X.  * j 1 * J £ ir ' _ x.  _ 1 ^ 1,1.  . i.  4 


with  experimental  results  of  Kumagai  et  al  that  show 
a reasonable  agreement.  Also  when  the  stoichiometric 
ratio  is  large,  it  is  found  that  the  droplet  lifetime 
is  shorter  than  the  flame  extinction  time,  and  an  ana- 
lytical expression  is  presented  giving  this  extinction 


time.  For  cases  in  which  the  extinction  time  is  much 


longer  than  the  lifetime  the  obtained  results  coincide 

A 
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with  those  obtained  by  Spalding  who  replaced  the  dro- 
plet by  an  instantaneous  source  of  vapor. 


1,2  Steady  State  Analysis  of  Hydrazine  Catalytic  Thrusters 


The  results  obtained  in  this  area  have  been  pre- 


sented in  detail  in  two  previous  reports  by  A.  Crespo, 
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E.  Fraga  and  A.  Munoz-'  and  A.  Crespo  . In  reference  5 


the  problem  is  stated  and  some  guidelines  for  the  solution 
are  given.  The  most  important  results  appear  in  ref. 6. 


The  model  examined  is  a onedimensional  stationary 
one  in  which  liquid  hydrazine  enters  the  reactor  and 
decomposes  yielding  enough  heat  to  vaporize  it  and  heat 
the  mixture.  The  hydrazine  decomposition  process  is 
assumed  to  follow  a two-step  reaction  mechanism  in  which 
hydrazine  decomposes  exothermically  yielding  ammonia, 
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hydrogen  and  nitrogen,  and  the  ammonia  in  turn  decomposes 
endothermically  producing  hydrogen  and  nitrogen.  The 
purpose  of  the  work  is  to  find  simple  expressions  for 
the  temperature  and  mass  fraction  distributions  along 
the  chamber  and  in  particular  at  the  exit  as  a function 
of  the  operating  conditions.  From  these  exit  conditions 
the  specific  impulse  that  would  be  obtained  if  the  mixture 
would  be  expanded  to  the  vacuum  is  obtained. 

Extensive  numerical  solutions  of  the  problem  are 
given  by  Kesten^  in  which  a large  number  of  parameters 
should  be  used  as  input  for  his  computer  programm  although 
he  also  gives  a numerical  correlation  of  his  results. 

g 

In  a work  carried  out  previously  at  INTA  also  numerical 
solutions  were  obtained  but  the  number  of  parameters 
was  practically  reduced  to  one,  although  the  results 
were  only  appliable  for  Shell  405  Catalyst  with  established 
values  of  the  velocity  and  stoichiometry  of  the  reaction. 

In  the  work  carried  out  with  the  present  Grant 
some  simplifications  were  introduced  to  the  model  of 
ref. 8 the  most  important  one  being  to  consider  that 
the  ammonia  decomposition  is  only  important  at  temper- 
atures close  to  the  maximum  one  that  is  reached  in  the 
reactor.  By  means  of  this  simplification  an  analytical 
solution  is  obtained  that  gives  conditions  at  the  chamber 
exit,  in  particular  specific  impulse'for  mixture  expanded 
to  the  vacuum  and  fractional  ammonia  decomposition, 
as  functions  of  operating  conditions  and  chamber  charac- 
teristics for  different  types  of  catalyst  particles 
and  chamical  kinetics.  These  results  show  a reasonable 
agreement  with  those  of  Kesten  , and  as  a matter  of 
fact  our  analytical  expression  for  fractional  ammonia 
decomposition  reproduces  some  aspects  of  his  correlation 
of  numerical  results. 


1 . 3 Separate  influence  of  catalyst  and  liquid  hydrazine 
temperatures  on  vacuum  start  characteristics  of  a 
hydrazine  catalytic  reactor. 


During  the  last  part  of  the  working  period  cor- 
responding to  this  grant,  work  has  been  done  on  studying 
the  temperature  distribution  within  a catalyst  particle 
during  the  period  of  liquid  hydrazine  penetration  in- 
side this  particle,  when  the  initial  temperatures  cf 
hydrazine  and  catalyst  are  different. 

One  of  the  most  important  problem  areas  in  the 
field  of  hydrazine  catalytic  decomposition  for  satellite 
microthrusters , is  the  catalyst  attrition  originated 
by  the  large  overpressures  created  inside  the  particles 
during  long  ignition  delays.  This  ignition  delay  time 
is  strongly  dependent  on  the  initial  temperature  at 
which  this  ignition  takes  place.  It  is  the  purpose  of 
our  work  to  study  the  separate  influence  of  both  the 
initial  temperatures  of  catalyst  and  liquid  hydrazine 
on  the  ignition  delay. 
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towards  the  particle  surface.  A different  model  has 
been  proposed  by  Webber  and  Schechter^  in  which  there 
is  not  definite  liquid  front,  but  gaseous  hydrazine 
and  products  are  allowed  to  coexist  with  liquid  hydrazine 
throughout  the  catalyst  particle;  although  this  may 
be  a more  realistic  model  than  the  one  used  here,  we 
shall  not  use  it  since  it  is  thought  that  a physical 
insight  of  the  process  can  be  obtained  without  recurring 
to  this  more  complicated  model. 

In  all  the  abovementioned  works  the  calculations 
were  performed  for  equal  or  nearly  equal  initial  tem- 
peratures of  liquid  hydrazine  and  catalyst. 

The  calculations  carried  out  in  this  report  show 
the  existence  of  a thermal  front  moving  behind  the  liquid 
front.  Ahead  of  the  thermal  front  the  liquid  hydrazine 
is  at  the  initial  temperature  of  the  catalyst,  Tq,  and 
behind  at  the  initial  temperature  of  the  liquid, 

{figure  2).  The  internal  structure  of  this  thermal  front 
has  been  calculated.  The  reason  for  the  existence  of 
the  front  lies  in  the  fact  that  convective  heat  transfer 
is  more  important  than  heat  conduction;  this  has  also 
been  observed  by  Webber  and  Schechter  Since  the  evap- 
oration and  reaction  of  hydrazine  takes  place  at  the 
liquid  front,  in  the  region  inside  it,  the  results  of 
the  report  would  indicate  that  the  only  temperature 
that  • would  have  an  influence  on  the  ignition  delay 
would  be  the  initial  particle  temperature  and  not  the 
liquid  hydrazine  temperature.  Some  experiments  carried 
out  by  Fouche  et  al^  confirm  that  the  ignition  time 
is  independent  of  liquid  hydrazine  temperature  if  the 
catalyst  bed  temperature  is  above  a certain  critical 
value,  however  for  smaller  values  of  the  catalyst  tern- 
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perature  the  ignition  time  becomes  dependent  on  both 
liquid  and  catalyst  temperatures,  thus  invalidating 
our  theory  in  this  second  situation.  Our  interpretation 
is  that  for  low  temperatures  the  ignition  time,  #.nd 
hence  the  residence  time  of  the  liquid  inside  the  par- 
ticle, becomes  long,  thus  decreasing  the  importance 
of  the  convective  effects  that  become  of  the  same  order 
of  the  heat  conduction.  The  limiting  conditions  under 
which  this  occurs  will  also  be  investigated. 

Other  experiments  on  this  area  have  been  carried 

12  13 

out  by  Sangiovanni  and  Kesten  and  Martney  and  Kesten, 
however  they  are  difficult  to  compare  with  our  theory 
because  they  either  use  the  same  temperature  for  liquid 
hydrazine  and  catalyst  or  do  not  state  clearly  the  dif- 
ference between  both  of  them. 

We  shall  liow  describe  in  more  detail  this  problem 
in  the  following  pages.  \ 

i 

i 

j 
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2.  TEMPERATURE  DISTRIBUTION  WITHIN  A CATALYST  PARTICLE 
DURING  PENETRATION  OF  LIQUID  HYDRAZINE 


The  problem  to  be  studied  consists  in  the  calculation 
of  the  temperature  distribution  within  a porous  spherical 
particle  that  is  initially  at  a temperature  Tq  when 
a liquid  at  a temperature  T^  penetrates  in  it  a prescribed 
flow  rate  Q that  is  a function  of  time.  The  main  assump- 
tions made  in  the  following  are  as  follows: 


i 


a)  Radial  symmetry 

b)  Liquid  hydrazine  penetrates  in  the  particle 
at  a prescribed  rate  Q(t)  that  we  suppose  to 

be  known.  We  shall  explain  later  how  this  datum 
can  br',  obtained  from  previous  calculations 
made  assuming  that  T^=Tq  > (refs.  8,9) 

c)  All  the  pores  between  the  liquid  front  radius  R 
and  the  particle  radius*  a, are  filled  with  liquid. 

d)  We  neglect  the  vaporization  rate  of  hydrazine 
compared  with  it.?  mass  flow  rate  during  this 
initial  unsteady  process.  This  may  be  justified 
by  comparing  typical  values  of  the  masses  of 
gas  and  liquid  ir.  the  particle  or  of  their 
densities  if  the  volumes  they  occupy  are  compa- 
rable; at  a typical  gas  pressure  of  100  atm. 
the  gas  density  is  an  order  of  magnitude  smaller 
than  t’  liquid  density.  Under  this  assumption 

Q and  R will  be  related  by: 


4 


t 

Tr(a3-R3)=  J Q(t)dt 
0 


i 
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(2) 


:^V^- 


Y being  the  volume  fraction  occupied  by  the  pores. 

e)  At  a given  time  and  location  within  the  parti- 
cle the  fluid  and  solid  temperature  are  equal. 

f)  There  is  an  overall  thermal  conductivity,  , 
that  is  supposed  to  be  temperature  independent. 

g)  Temperature  changes  due  to  other  causes  besides 
the  unequal  initial  temperatures  of  liquid 

and  particle  are  neglected.  Sangiovanni  ,'nd 
Kesten  found  in  their  numerical  calculations 
temperature  rises  of  less  than  4°K  for  cases 
in  which  these  initial  temperatures  were  equal. 

With  these  assumptions  the  equation  governing 
heat  transfer  within  the  particle  in  the  liquid  region, 
R<r<a,  is 


QplCl  3 T 
4n  r2""  3r 


— ~ — (Kp  r2~~— ) 
2 dr  p 9r 


pc  is  an  average  value  of  specific  heat  per  unit  volume 

pcts  Psds(1~Y)+P1c1T  <4) 

p and  p-  are  the  densities  of  the  solid  forming  the 
catalyst  and  of  the  liquid  hydrazine,  and  cg  and  c^ 

the  corresponding  specific  heats  per  unit  mass. 


2.1  First  Approximation 


It  is  obtained  by  neglecting  the  las  term  of  equa- 
tion (3)  representing  heat  conduction.  3y  comparing 
this  term  with  the  second  one  representing  heat  convec- 
tion, it  is  obtained  that  their  ratio  is  of  order  of 
the  ratio  of  the  liquid  residence  time  within  the  par- 
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tide,  tr>  and  the  characteristic  time  for  heating  of 
the  particle  by  conduction,  t 

4/3* a 3 plCla2 

r Q0  c“  Kp  (5> 

Where  Qq  is  a typical  value  of  q.  For  particles  of  radius 
of  the  order  of  .5mm.,  refs.  8,9  give  tor  t values  below 
0.1  sec.  for  temperatures  above  290°K.  A typical  value 
of  t for  that  particle  radius  is  1 sec. 

Equation  (3)  without  the  conduction  term  becomes 
of  the  hyperbolic  type  and  has  a general  solution  of 
the  form 

T-f(?)  (6) 

Where  £=  *(a3-R3) |-Tr(a3-r3)  (7) 

J YP  j^c  j 3 

and  equation  (2)  has  been  used. 

The  initial  and  boundary  conditions  to  be  satisfied 
are 


O 

II 

u 

r<a,  T=T  o 

(8) 

t>0 , 

r = a,  TaT'^ 

(9) 

With 

these  conditions  f(?) 

becomes  a step  function  such 

that 

HO- 

=T.  for  £>0 

(10a 

for  ^<0 


f (5)“T0 


(10b) 


This  means  that  there  is  a thermal  front  of  radius  R 

« 

given  by  £=0 


— j-tt  (a3-R3)  = pc  -|—Tt(a3-R3) 


behind  which  the  liquid  and  particle  temperatures  are 
equal  to  the  injected  hydrazine  temperature  T^,  and 
ahead  cf  it  equal  *;o  the  initial  particle  temperature 
T o « The  thermal  front  goes  always  behind  the  liquid 
front,  Rt>R,  as  can  be  seen  from  equation  (11)  (figure  i). 
Equation  (11)  expresses  a heat  balance  between  the  heat 
convected  by  the  liquid  hydrazine  inside  the  particle 
and  the  non-steady  particle  heating. 


2.2  Thermal  front  structure 


To  study  the  flame  front  structure  we  shall 
assume  it  to  be  tbin  compared  wi  h che  particle  radius  a, 
and  make  the  following  change  of  variables 


r=Rt(t)+6x 


where  6 is  a characteristic  thickness  of  the  front 


6<<a 


and  x a non-dimensional  length  of  order  unity.  \ ie 
time  is  made  nondimens ional  with  a quantity  o !'  the  order 
of  t„ 


l)o 

2UQJ  pc 


■piJSPBJS  V&L grass’ 


-13 
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By  making  this  change  of  variables  in  equation  (3), 
it  is  obtained 


3T  , of  . 9T 

sf  + P(T)x  u 


3ZT 

3x2 


'(15) 


after  neglecting  terms  of  order  6/a,  and  choosing  6 
of  the  form 


$2-  _^E  ..Zira? 

P j c2  Qo 


(16) 


and 


_5_ 

Qo 


07) 


and  using  equations  (11)  and  (2)  for  the  thermal  front 

radius  equation  (16)  expresses  that  condition  (13) 

is  satisfied  if  t is  much  smaller  than  t (as  given 

r c 

by  equation  (5)  ) . 


The  appropiate  boundary  and  initial  conditions 
for  equation  (15),  if  condition  (13)  is  satisfied,  are 


X-*co 

r=T0 

(18a) 

x->-» 

T=Tj 

U'ib) 

t = 0 

x<0  T=T  o 

(19) 

Equation  (15)  and  conditions  (18)  and  (19)  admit  a 
similarity  solution  of  the  form 


I 


1 


1 


3^SS**l5?M'WWS,~yj.  Jr^.45^^,j^.jyjyi.'ijjjj^  ***  ••  r^T^rrsr^y: 
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T-TJ+ 


T0“T. 


erf  c z 


(20) 


Where 


z= 


$ (t) 


(21) 


and  (j>(t)  is  a function  of  time  satisfying  a differential 
equation  of  the  form 


^-3$2-l=0 


(22) 


with  the  initial  condition 


tb0  $=0 


(23) 


and  3 a funccion  of  t given  by 


Qo 


R. 


(24) 


The  solution  of  (22)  and  (23)  is  gi-^en  by 


T T T 

<j>2  = 2exp(  J 23dx  ) J exp(-J  23dx)dx 
0 0 0 


(25) 


and  using  the  expression  (17)  for  3 and  equations 
(2)  and  (11)  for  Rt>  it  is  obtained  for  <{> 


x R H 

/ L.  v 


*2=  2<-f->4 


0 a 


(26) 


The  thermal  front  has  a thickness  that  increases  with 
time  according  to  '•he  law  given  ’jy  equation  (26) 


•M^jML&iuiii^iases^'^^^uuM 
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2.3  Conclusions  and  Discussion 


The  present  theory  indicates  that  when  there  is 
a situation  for  which  condition  (13)  is  sacisfied  the 
appropiate  temperature  to  take  into  consideration  to 
calculate  the  liquid  residence  time  within  the  catalyst 
particle  is  the  initial  particle  temperature,  aud  the 
liquid  hydrazine  temperature  does  not  have  influence 
on  this  starting  period.  As  the  liquid  residence  time 
increases  the  thickness  of  the  thermal  front  increases 
and  may  become  of  the  order  of  the  particle  radius. 

This  limiting  situation  arises  when 


6 -a  , 6ca2  . 

~KP  *r  - 

An  analytical  expression  for  tr 
which  it  is  shown  that  t behav 
and  temperature  in  the  form 


(27) 

is  given  in  ref.  8 in 
es  with  particle  radius 


a 

— exp 


, 5809. 

\ fp  ) 9 


T in  °K 


(28) 


showing  that  thr’s  criteria  is  particulary  sensitive  to 

temperature  variations.  For  initial  particle  temperatures 

below  290°K  and  particle  radius  of  0.5mm.  t becomes 

r 

also  dependent  on  liquid  hydrazine  temperature,  according 
to  the  results  of  refs. 8,9,  that  give  then  for  t values 
above  -lsec. 


A different 
what  happens  for  1 
smaller  of  the  ord 


situation 

to 

be 

discussed 

is  also 

arger  valu 

es 

of 

T when  t 

becomes 

er  of  the 

mills 

econd  for 

a=0.  5mm. 

> 
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then  the  results  of  refs.  8,9  ahow  that  the  liquid  penetrates 
a small  depth,  "e”,  inside  the  particle  and  criteria  (27) 
should  be  substituted  by 


See1* 


K t 
P r 


(29) 


However,  an  asymptotic  analysis  carried  out  in  ref. 8 

-1/2 

shows  that  "e"  decreases  precisely  like  t , so  that 

both  effects  compensate  and  the  above  presented  theory 
is  expected  io  be  valid  in  this  limit. 
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